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Abstract: As Lindzen (1981) had shown, small-scale gravity waves (GW) produce the observed 
reversals of the zonal-mean circulation and temperature variations in the upper mesosphere. The 
waves also play a major role in modulating and amplifying the diurnal tides (DT) (e.g., Waltershcid, 
1981; Fritts and Vincent, 1987; Fritts, 1995a). We summarize here the modeling studies with the 
mechanistic numerical spectral model (NSM) with Doppler spread parameterization for GW (Hines, 
1997a, b), which describes in the middle atmosphere: (a) migrating and non-migrating DT, (b) 
planetary waves (PW), and (c) global-scale inertio gravity waves. Numerical experiments are 
discussed that illuminate the influence of GW filtering and nonlinear interactions between DT, PW, 
and zonal mean variations. 

Keywords: Theoretical modeling. Middle atmosphere dynamics. Gravity wave interactions, Migrating 
and non-migrating tides, Planetary waves. Global-scale inertio gravity waves. 

1. Introduction 

The fundamental properties of diurnal tides (DT) and planetary waves (PW) are well understood 
(Chapman and Lindzen, 1970; Holton, 1979; Kato, 1980; Volland, 1988). Measurements from the 
ground and with spacecraft have shown that the DT in the mesosphere and lower thermosphere exhibit 
large seasonal and inter-annual variations, and that they are modulated by PW (e.g., Avery et al., 1989; 
Manson et al., 1989; Vincent et al., 1989; Gille et al., 1991; Hays et al., 1994; Burrage et al., 1995a, b; 
McLandress et al., 1996; Wu et al., 1998; Leblanc et al., 1999a, b; Huang and Reber, 2003, 2004; Wu 
and Jiang, 2005; Zhang et al., 2006; Muktarov et al., 2009; Huang et al., 2006, 2010). 

In the mesosphere, the DT, PW, and small-scale gravity waves (GW) attain large amplitudes, and 
their interactions have been the subject of theoretical studies (e.g., Walterscheid, 1981; Fritts, 1984, 
1995a, b; Fritts and Vincent, 1987; Vial and Forbes, 1989; Forbes et al., 1991; Forbes, 1995). 
Numerical models have been employed to study the DT under the influence of the zonal-mean 
circulation, eddy viscosity, and GW interactions (e.g., Lindzen and Hong, 1974; Forbes and Hagan, 
1988; Forbes and Vial, 1989; Akmaev et al., 1996; Miyahara and Forbes, 1991; Hagan, 1996; Hagan et 
al., 1993, 1995, 1999a; Yudin et al„ 1997; McLandress, 1997a, 2002a, b; Miyahara et al., 1993, 1999; 
Geller et al., 1997; Mayr et al., 1998; Akmaev, 2001a, b). 

Mayr et al. (2010), referred to as MMCH, summarize the zonal-mean variations in the middle 
atmosphere, which are generated in the numerical spectral model (NSM) with parameterized GW. In 
this companion paper, we discuss the tides, planetary waves, and global-scale inertio gravity waves. 

Chan et al. (1994, 1995) introduced the NSM, which is fully nonlinear and is integrated from the 
surface into the thermosphere. Formulated in terms of vector spherical harmonics, the model allows us 
to study the individual and interacting dynamical components with zonal wave numbers m = 0 to 4, as 
illustrated in Fig. 1 of MMCH. The zonal-mean (m = 0) heating rates for the middle atmosphere are 
taken from Strobel (1978), and the excitation rates for the migrating DT in the troposphere and 
stratosphere are taken from Forbes and Garrett (1978). Solar EUV heating is applied in the 
thermosphere. The planetary waves (PW) are generated internally by instabilities. Newtonian cooling 
describes radiative loss (Wehrbein and Leovy, 1982; Zhu, 1989). Starting with Mengel et al. (1995), 
the NSM has been ran with the Doppler spread parameterization (DSP) for GW developed by Hines 
(1997a, b). In the DSP formulation, GW momentum is strictly conserved, which is important for wave 
filtering and the associated nonlinear interactions. The DSP comes with uncertainties that affect the 
height dependent GW momentum source and eddy viscosity, and the input parameters were changed 
over the years as documented in the literature. In all DSP applications, for simplicity, the GW source 
at the initial height was chosen to be isotropic, time independent, and hemisphericatly symmetric. 

By way of introduction, wc present in Fig. 1 numerical results from 2D (a) and 3D (b) versions of 
the NSM. Applying the same GW 7 source at the initial height, the zonal-mean (m = 0) winds are 
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presented, which describe the quasi-biennial oscillation (QBO) and semi-annual oscillation (SAO) near 
the equator. From the 2D model (Fig. 1 a), the wind velocities are much larger than those in 3D (Fig. 
lb). The GW in 2D are solely available to generate and amplify the QBO and SAO. In the 3D model, 
the GW apparently amplify the tides and planetary waves at the expense of the equatorial oscillations. 
This illustrates how GW interactions can influence, and couple, the dynamical components of the 
atmosphere. 

In the following we show how GW processes affect the tides and planetary waves as enumerated in 
Fig. 2, and we present numerical experiments that provide understanding. 

2. Diurnal Tides 

Generated primarily by solar heating, the dominant diurnal tides propagate with the Sun westward 
and are referred to as migrating tides. The non-migrating tides, not propagating with the Sun, are 
generated in the NSM with planetary waves (PW) — and they are later discussed. Here we deal 
primarily with migrating tides, but the oscillations include also non-migrating components. 

The basic effect of GW momentum deposition on the fundamental 24-hour diurnal tide is 
demonstrated in numerical results taken from Mayr et al. (1998). Fig. 3 shows the amplitude and 
phase of the meridional winds, generated with zonal wave number m = 1 . Computer solutions are 
presented with (a) and without (b) GW source. The differences are small below about 60 km. But at 
higher altitudes, the GW increase the amplitudes of the tide and reduce its vertical wavelengths — as 
illustrated in Fig. 4. 

With Hines' GW parameterization, the waves amplify the tides, as shown by McLandress (1997b). 
Due to the nonlinear positive feedback associated with critical level absorption, the GW deposit 
momentum in vertical wind shears (e.g., Walterscheid, 1981). This is evident in Fig. 6 of MMCH, 
which describes the momentum source that generates the QBO. The kind of wave interaction required 
to produce the QBO is amplifying the diurnal tide. 

2.1. Seasonal Variations with Planetary Wave Modulation 

Satellite measurements show that the diurnal tide (DT) in the upper mesosphere exhibits large 
seasonal variations with peak amplitudes near equinox (Hays et al., 1994; Burrage et al., 1995a, b; 
McLandress et al., 1996; Huang and Reber, 2003). The observations are reproduced in the NSM 
(Mayr et al., 1998, 1999), which describes the DT and the zonal-mean variations, including QBO and 
SAO, generated with GW interactions. This is shown in Fig. 5, where we present numerical results for 
the zonal-mean winds (a) and the meridional wind amplitude of the tide (b). The DT in Fig. 5b is 
shown at 95 km and 18° latitude where the meridional winds peak. The mean zonal winds in Fig. 5a at 
50 km are eastward in winter and westward in summer due to solar heating, and at 80 km the winds 
reverse direction due to GW filtering, as discussed in MMCH. 

Our model results in Fig. 5 show that the peak amplitudes of the tide occur where the zonal winds 
change direction, which indicates that GW filtering is involved. The upward propagating GW interact 
with the large zonal winds during summer and winter months, and deposit momentum. Less GW 
momentum is thus available to amplify the tide around solstice (January/July), which could produce 
the amplitude maxima during equinox that are observed. 

As Fig. 2 illustrates, the dynamical situation is more complicated than suggested and portrayed in 
Fig. 5. Numerous studies have shown (e.g., Lindzen and Hong, 1974; Forbes and Hagan, 1988; Hagan 
et al., 1993) that the zonal mean temperature and wind fields affect the diurnal tides and their seasonal 
variations. Moreover, the seasonal variations of the planetary waves (PW) absorb and in turn modulate 
the GW, which contributes to the variability of the tides. GW filtering by the diurnal and semi-diurnal 
tides can produce interactions that affect their seasonal variations. Mayr et al. (2001) discussed 
numerical experiments that describe these processes, and in the following we present some model 
simulations from that study. 
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The numerical results in Fig. 6 show the amplitude and phase of the diurnal tide (DT) at 95 km 
and 18° latitude. The standard solution (solid line, large + marks) accounts for all dynamical 
interactions between tides, PW, and zonal mean variations under the influence of GW. For 
comparison, a numerical experiment is presented (dotted lines, small + marks), in which the GW are 
turned off, but only for the tides and PW (m > 0). In this case, the GW still influence the DT through 
the zonal-mean wind and temperature variations — and that interaction produces seasonal variations 
with equinoctial maxima (dotted lines). With GW source (solid lines), however, the equinoctial 
maxima are much more pronounced, and the phase is shifted to earlier local times. The interaction of 
the DT with the zonal mean circulation contributes to the semi-annual seasonal variations, shown by 
McLandress (1997a), but GW filtering is more important. 

As shown in Section 5, the planetary waves (PW) are generated in the NSM internally by 
baroclinic instabilities, which can develop in the zonal-mean temperature variations, as proposed by 
Plumb (1983) and Pfistcr (1985) and confirmed in modeling studies (e.g., Norton and Thubum, 1996, 
1999). Fig. 7 shows the PW for m = 3 at 80 km, computed with and without GW interaction. The PW 
appear primarily around summer/winter solstice where instabilities develop. And the waves are greatly 
amplified by GW, which explains the enhanced PW modulation of the DT in Fig. 6. 

The PW also affect the seasonal variation of the DT. This is shown in Fig. 7b, where we present 
the 5-day mean of the meridional wind amplitudes for the tide at 95 km, computed with and without 
the PW for m = 3 and 4. Without PW, the DT is much larger during summer months. This indicates 
that GW filtering by the PW is involved. The GW momentum expended to amplify the PW is not 
available to amplify the DT at higher altitudes. 

The semi-diurnal tide (SDT) peaks at mid latitudes, and in Fig. 8a we present the computed zonal 
wind amplitudes at 95 km that are produced with and without GW. With the GW source, the PW 
modulations are large and similar to those of the DT in Fig. 6. The amplitude of the SDT is amplified 
by GW, and it peaks during winter and summer months when the DT is weakest. This relationship 
between the tides is made evident in Fig. 8b, where we present the 5-day mean of the SDT computed 
with and without the excitation of the DT. Without the DT, the amplitude of the SDT tends to be 
larger, and the differences are largest during the spring and fall seasons. Around equinox, the GW 
apparently amplify the DT at the expense of the SDT. (Without GW forcing, Mayr et al., 2001 report 
that the SDT is not much affected by the absence of the DT). 

2.2. QBO-generated Inter-annual Variations 

Satellite and ground-based measurements show that the diurnal tide (DT) reveals relatively large 
inter-annual variations in the upper mesosphere (e.g., Burrage et al., 1995a; Lieberman, 1997; Vincent 
et al., 1998). Hagan et al. (1999a) carried out an analysis with the Global Scale Wave Model (GSWM) 
and showed that the seasonal variations of the DT are significantly affected by linearized advection 
associated with the mean (m = 0) zonal winds. Advection produced with the observed QBO, however, 
could not reproduce the inter-annual variations of the DT (Hagan et al., 1 999b). Applying Hines GW 
parameterization, McLandress (2002c) generated an 18-month QBO with zonal winds close to 40 m/s 
at 30 km. His analysis demonstrated that significant inter-annual variations were produced by the 
advection term from the meridional gradient of the QBO mean zonal winds. 

The numerical results presented here are taken from Mayr and Mengel (2005). Fig. 9 shows the 
computed stratospheric QBO, and the inter-annual variations of the DT in the upper mesosphere. Near 
the equator, the zonal-mean (m = 0) zonal winds at 30 km (Fig. 9a) vary with a period close to 24 
months, and that periodicity is evident in the amplitude of the DT at 95 km (Figs. 9b, c). A 
pronounced signature of inter-annual variations is also seen in the DT at 18° latitude (Fig. 9d), where 
amplitude maxima develop during the spring/fall seasons as observed and discussed in Section 2 . 1 . 

The relationship between the stratospheric QBO and inter-annual variation of the DT is clearly evident, 
and the question is what are the mechanisms that produce this connection. 
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In agreement with Hagan et al. (1999a, b) and McLandress (2002c), our analysis with the NSM 
confirms the importance of advection for the inter-annual variation of the diurnal tide. But it also 
shows that the QBO-modulated GW momentum source dominates in the upper mesosphere. 

The NSM is nonlinear, and we present here the linearized advection terms for the meridional and 
vertical components of the DT, 


, vSU , 

(v) — +(w) — 
w r<98 w Sr 


( 1 ) 


where U is the mean (m = 0) QBO zonal wind; (v), (w) represent the time-average amplitude of the 


meridional and vertical winds of the DT, respectively; and r, 0 are Earth radius and latitude. (The 
zonal advection term can be ignored.) 

The numerical results for advection are shown in Fig. 10. For the mean zonal winds, a 1-year 
r unnin g mean was taken to describe the QBO variability, and the terms were normalized to represent 
% variations in relation to the acceleration of the tide. Since meridional advection (Fig. 10a) becomes 
important away from the equator, the term is presented at 1 8° latitude where the DT peaks. The 
vertical advection (Fig. 10b) is shown at 4° near the equator, where the QBO zonal winds peak. 
Comparable in magnitude, the meridional and vertical advection terms reveal pronounced QBO like 
modulations around 30% in the stratosphere. But in the mesosphere above 60 km, the 24-month QBO 
signatures almost disappear. With QBO zonal winds of 20 m/s in the lower stratosphere, the advection 
tenns cannot reproduce the inter-annual variations of the tide that are generated with the NSM. 

In the upper mesosphere, GW become important. This is shown in Fig. 11, where we present the 
QBO zonal winds (a), the GW momentum source (b), and the variability of the tide at 95 km (c). For 
the amplitudes of the GW source (Fig. 1 lb), the normalized relative variations (%) are presented, 
analogous to Fig. 10. Above 80 km, year-to-year variations of about 20% are generated, and the 
period is close to 24 months of the QBO. The yearly variation of the DT (Fig. 11c) shows amplitude 
maxima and minima that are aligned with the GW source (Fig. 1 lb), which demonstrates that the GW 
modulate the tide. The GW are apparently filtered by the QBO zonal winds in the middle atmosphere. 
During the years 3, 5, 7, the waves appear to be preferentially absorbed, and the resulting momentum 
deficits reduce the amplitudes of the DT at higher altitudes. 

In the stratosphere, the momentum advection from the QBO zonal winds is very effective in 
modulating the tide. But in the mesosphere, the GW filtered by the QBO attain large amplitudes and 
provide the momentum source that generates the inter-annual variations. 


3. Non-migrating Tides 

The so-called diurnal and semi-diurnal migrating tides (MT) propagate westward with the Sun, 
occupy m = 1 and 2, and have periods of 24 and 12 hours, respectively. The non-migrating tides 
(NMT) also have periods of 24 and 12 hours, commensurate with solar heating. But unlike the MT, 
the NMT can propagate eastward, be stationary (m = 0), or occupy any wave number. 

NMT have been observed in the middle atmosphere with satellite borne wind and temperature 
measurements (e.g., Lieberman, 1991; Forbes et ah, 1995, 2003; Talaat and Lieberman, 1999; 
Oberheid et al., 2002; Huang and Reber, 2004). The observed phase structure of the NMT indicates 
that the tides propagate upward (e.g., Lieberman. 1991; Talaat and Lieberman, 1999). 

Hagan and Forbes (2003) demonstrated that NMT can be generated by topography and latent heat 
release associated with convection. They showed that this tropospheric source contributes 
significantly to the observed tidal amplitudes in the upper mesosphere. 

Teitelbaum and Vial (1991) proposed that nonlinear coupling between migrating tides (MT) and 
planetary waves (PW) can generate NMT. Observed correlations between MT and PW activity 
provide observational evidence for this mechanism (e.g., Forbes et al., 1995; Talaat and Lieberman, 
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1 999, Pancheva et al., 2002); and modeling studies have shown that the interaction can produce NMT 
with large amplitudes (e.g., Miyahara et al., 1997, 1999; Hagan and Roble, 2001, Coll and Forbes, 
2002 ). 

The NSM generates the NMT with PW (Mayr et al., 3003, 2005a, b), and Fig. 12 presents the 
meridional wind amplitudes for the 24-hour tides at 100 km. As discussed in Section 2.1, the dominant 
diurnal westward propagating MT with m = 1 peaks near 20° latitude (Fig. 12a). The amplitudes are 
largest (~ 100 m/s) during the spring and fall seasons, owing in part to GW filtering by the mean zonal 
winds. The m = 1 eastward propagating NMT (Fig. 13b) is highly variable, and it peaks around the 
equator with velocities between 10 and 30 m/s. Similar amplitudes are also generated in the stationary 
NMT (m = 0) shown in Fig. 12c. 

The mechanism for generating NMT with PW is illustrated in Fig. 13. Numerical results are 
presented for the stationary m = 0 NMT, which is produced through non-linear interactions between 
PW and the diurnal MT, each having wave number m = 1. Fig. 13a represents the standard solution for 
the meridional wind amplitude of the NMT and the zonal wind PW at 1 00 km and 1 8° latitude. During 
the initial 6 months of the model ran, the amplitudes are relatively small. But later, the PW and NMT 
attain values that approach 20 and 40 m/s respectively. The amplitude variations show evidence of 
correlations. For comparison, we present in Fig. 13b numerical results from a model run without the 
zonal-mean heat source that produces the PW through instabilities. After the initial start-up noise, the 
PW and NMT are no longer generated. 

The NSM is fully nonlinear, and a number of processes could produce the interaction between PW 
and MT that generates the NMT (Figs. 12 and 13). Gravity wave (GW) filtering can also produce 
nonlinear interactions, as shown for example in MMCH for the 5-year oscillation, and in Sections 2.2 
for the inter-annual variations of the diurnal tide. For the NMT, the proposed nonlinear GW 
interaction is illustrated in Fig. 14. Upward propagating GW encounter a long-period PW and deposit 
eastward momentum in the eastward phase of the wind oscillation. The wave momentum given to the 
PW is then missing in the GW when they encounter the MT at higher altitudes. Since the GW amplify 
the tide, its eastward (positive) phase is weaker than the westward phase. And GW filtering by the 
negative phase of the PW produces the weaker negative phase of the tide. In effect, the tide is 
modulated by PW, which is a signature of nonlinear interaction. 

Fig. 15 is taken from Mayr et al. (2005b) and presents numerical results for the m = 0 NMT from 
solutions with and without GW source, which indicates that the waves contribute significantly to the 
nonlinear coupling between MT and PW. In this study, the GW source was simply turned off after 5 
years, and the subsequent 3 months are shown for the meridional and zonal winds at 1 00 km. Having 
integrated the model sufficiently long, the idea was that the inertia in the zonal-mean temperature and 
wind velocity could maintain, for a limited time, the PW and MT at full strengths. For the PW, our 
numerical results confirm that. But for the MT involved, the zonal and meridional winds vary in 
opposite directions when the GW source is turned off. The large differences apparent in Fig. 15 thus 
must be taken with some caution. 

A more definitive demonstration of the GW process is presented in Fig. 16, taken from Mayr et al. 
(2003, 2005b). The focus here is the nonlinear GW interaction specifically, which is involved in 
generating the m = 0 NMT. The model was ran for perpetual equinox, without time dependent solar 
heating. For comparison, the model was then integrated again but with a GW momentum source for m 
= 0 that was forced to remain constant after 1 8 months. In this control experiment, only the non-linear 
GW interaction is suppressed that can generate the NMT. All the other GW processes are retained that 
filter the MT and produce the PW through instabilities. The building blocks of the NMT are generated 
without interference. 

Fig. 1 6a shows the computed meridional wind amplitude of the westward propagating m = 1 
diurnal tide at 100 km. Both solutions are presented, and the differences are small. For the m = 1 PW 
in Fig. 16b, the time independent GW source produces somewhat larger velocities. Fig. 16c then 
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shows the m = 0 non-migrating diurnal tide computed with and without the time independent GW 
momentum source, and the differences are large after 18 months. Without the GW interaction (lower 
panel), the NMT is greatly reduced, although the PW involved are somewhat stronger (Fig. 16b). This 
clearly demonstrates that the GW contribute significantly to the nonlinear coupling between MT and 
PW. 

4. Planetary-scale Inertio Gravity Waves 

In the upper mesosphere, oscillations with periods around 1 0 hours have been observed with 
different measurement techniques. Analyzing radar wind measurement at Saskatoon (52° N), Manson 
et al. (1982) and Manson and Meek (1990) discovered persistent oscillations with periods of 9-10 and 
16 hours, which were attributed to nonlinear coupling between the semidiurnal tide and planetary 
waves. Similar kind of variations were also reported by Ruster ( 1 994) and Sivjee et al. ( 1 994). In 
optical measurements of winds and temperatures at 90 km, Hemadez et al. (1992, 1993) observed 
pronounced westward propagating waves with a period of 10. 1 hours near the South Pole. The 
inferred large vertical wavelength indicated that a Lamb wave was observed. At the South Pole (88° 
S), westward propagating wind oscillations with periods 7-12 hours were observed with radar 
measurements (Forbes et al., 1999; Portnyagin et al., 2000). Lamb waves and "normal modes" were 
invoked, but inertio gravity waves, nonlinear coupling between tides and planetary waves were also 
suggested as possible explanations. With measurements from a number of instruments at northern 
high latitudes, Wu et al. (2002) observed 10-hour oscillations in the temperature and wind fields of the 
upper mesosphere and lower thermosphere. A zonal wave number m = 5 was inferred, and the vertical 
wavelength was estimated at 53 km. The authors proposed that the oscillations could be generated by 
the nonlinear interaction between the semidiurnal tide and 2-day planetary waves. 

The NSM generates in the upper mesosphere oscillations with periods around 10 hours having 
zonal wave numbers m = 0 - 2 (Mayr et al„ 2003, 2004). In Fig. 17, the computed meridional wind (a) 
and temperature (c) variations are shown at 84° N for a short time span in March of the third model 
year when the oscillations are relatively large. Like diurnal tides, the waves propagate upward, 
increase in amplitude, and the vertical wavelengths are around 20 km. The winds and temperature 
perturbations attain at 100 km values close to 20 m/s (Fig. 17a) and 5 K. (Fig. 17c), respectively. The 
power spectrum for the meridional winds (Fig. 17b), obtained from a time span of 20 days, shows that 
the amplitudes for the 9 and 10-hour waves exceed those of the non-migrating 12-hour semidiurnal 
tide. As was shown for the diurnal tides, the large wave amplitudes in the mesosphere are presumably 
produced in part by GW momentum deposition due to critical level absorption in vertical wind shears. 

In our model results, the temperature oscillations for m = 0 tend to peak at the poles. The m = 1 
horizontal winds can also peak at the poles, but that depends on the direction of propagation as shown 
in Fig. 18. Presented here are the meridional wind spectra at 100 km obtained from a 6-month time 
span following December solstice. With a running window of 1 5 days, the westward (red) and 
eastward (blue) propagating waves were derived. As expected, the westward propagating 12-hour 
semidiurnal tide dominates, which attains peak amplitudes close to 40 m/s near the poles in both 
hemispheres. The spectra also show pronounced amplitudes for oscillations with periods close to 10 
and 9 hours, which propagate westward and reach amplitudes exceeding 20 m/s near the North Pole. 
The weaker eastward propagating waves have amplitudes ~ 10 m/s, confined to mid latitudes. 

The dynamical conditions characterizing these oscillations have much in common with the earlier 
discussed diurnal tides. With periods close to 12 hours, the Coriolis force is important; the vertical 
wavelengths have similar magnitudes; and the oscillations have low zonal wave numbers (m < 4). 
Based on classical tidal theory (Longuett-Higgins, 1968), close examination shows that the oscillations 
are related to Class- 1 waves or Hough modes, which have also been referred to as Inertio Gravity 
Waves (Andrews et al., 1987: Volland, 1988). Since the waves occupy low zonal wave numbers, we 
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refer to them as Planetary-scale Inertio Gravity Waves (PIGW). The P1GW differ from typical Inertio 
Gravity Waves (1GW), which have also periods between 5 and 11 hours but are observed with much 
shorter horizontal wavelengths from 500 to 3000 km (Fritts, 1984; Fritts et ah, 1984; Muraoka et ah, 
1988; 1994; Tsuda etah, 1989). 

With relatively long horizontal wavelengths and short oscillation periods, the P1WG in our model 
have phase velocities in excess of the background zonal winds. Unlike the slower planetary waves 
(PW) with longer periods, or IGW with shorter wavelengths, the fast PIGW can propagate through the 
atmosphere without dissipation due to critical level absorption. 

The seasonal variations of the PIGW are shown in Fig. 1 9a for the m = 1 meridional winds at 100 
km. In both hemispheres, the waves are generated preferentially around winter solstice, extending into 
spring. Consistent with Fig. 18, the westward and eastward propagating oscillations are respectively 
confined to high and mid latitudes. Our analysis shows that the waves arc not generated when the 
zonal-mean heat source of the atmosphere is turned off. Like the PW, later discussed, the PIGW are in 
part generated by instabilities. 

The diurnal tides are also involved. And this is shown in Fig. 1 8b, where we present the m = 1 
PIGW computed without the excitation of the migrating diurnal tides. The amplitudes are of 
comparable magnitudes, but the seasonal variations differ considerably. The PIGW in Fig. 1 8a are 
apparently generated like the non-migrating tides through nonlinear coupling between PW and 
migrating tides, as had been suggested in several papers (e.g., Manson et al., 1982; Manson and Meek, 
1990; Forbes et ah, 1999; Portnyagin et ah, 2000). 

In light of the results shown in Fig. 1 8b, it is not surprising that 12-hour pseudo tides are generated 
without the solar excitation of migrating tides (Talaat and Mayr, 2010). In the upper mesosphere, and 
in the NSM, this class of non-migrating tides contributes significantly to the semi-diurnal oscillations 
at high latitudes. 

5. Planetary Waves 

The dynamical properties of planetary waves (PW) are well understood and have been discussed in 
the textbook literature (e.g., Holton 1979; Volland, 1988). Topographic forcing and tropical 
convection can produce the PW originating in the troposphere, and the baroclinic instability has been 
invoked for the PW generated in the mesosphere (Plumb, 1983; Pfister, 1 985). As the PW propagate 
through the middle atmosphere, they attain large amplitudes in the upper mesosphere as observed with 
ground-based and spacecraft measurements (e.g., Muller and Nelson, 1978; Craig and Elford, 1981; 
Burks and Leovy, 1986; Tsuda et al., 1988; Phillips, 1989; Poole, 1990; Wu et al., 1993; Fraser et al., 
1993; Clark et al., 1994; Fritts and Isler, 1994; Harris, 1994; Smith, 1996; 1997; Meek et al., 1996; 
Deng et al., 1 997; Fritts et al., 1999; Talaat and Lieberman, 1999; Wang et al., 2000; Talaat et al., 
2001 ). 

In the NSM, the PW are generated internally by instabilities associated with the zonal-mean 
variations of temperature, density, and zonal winds. As shown in Fig. 13, the PW are not excited when 
the zonal-mean heating is turned off. The model does not describe the PW that are produced in GCMs 
by topography and tropical convection (e.g., Giorgetta et al., 2002; Horinouchi et al., 2003). 

The zonal-mean (m = 0) temperature variations in the NSM are presented in Fig. 20a for December 
solstice, and Fig. 20b illustrates the associated dynamical processes. Owing to solar heating in the 
stratosphere and lower mesosphere, the temperatures decrease in winter and increase in summer to 
produce respectively eastward (positive) and westward (negative) zonal winds. Generated by GW 
interaction with the zonal circulation (Lindzen, 1981), the zonal winds reverse direction in the 
mesosphere, and in geostrophic balance the temperatures him colder in summer than in winter. Due to 
GW, the meridional circulation in the mesosphere increases to such magnitude that more energy is 
removed in summer than provided by solar radiation. 
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Around the tropopause, cumulus heating at low latitudes is generating zonal jets that are observed 
at mid latitudes. This is simulated in the NSM with a time independent tropospheric zonal-mean heat 
source. The heat source generates in the model the temperature variations in Fig. 20, which peak at the 
equator in the troposphere but increase to higher latitudes in the lower stratosphere, in qualitative 
agreement with observations. 

The dynamical conditions in the troposphere are to some extent similar to those occurring in 
middle atmosphere, as illustrated in Fig. 20b. In the mesosphere, GW filtering causes the zonal winds 
and temperature variations to reverse at higher altitudes, which is conducive for generating PW 
through baroclinic instabilities (Plumb, 1983). We believe that GW interactions could also produce the 
reversing temperature variations around the tropopause, to generate, in the NSM, the tropospheric PW 
through instabilities. 

In Fig. 21, we present at mid latitudes the PW for m = 1 and 3 at 90 km, which are generated with 
(a) and without (b) the tropospheric heat source. With the tropospheric source, the long-period PW 
amplitudes for m = 1 are highly variable but tend to peak during the spring and fall seasons. Without 
the source, the oscillations are greatly reduced and much less variable; and the differences for the m = 

2 PW show similar patterns. For the PW with m = 3, the tropospheric source also increases the wave 
amplitudes but mainly in winter and around equinox. These short-period waves are generated mainly 
in the stratosphere and mesosphere. 

The seasonal/latitudinal variations of the m = 3 PW are shown in Fig. 22 at 60 and 90 km. The 
amplitudes are computed with running windows, in which the oscillation periods are multiplied by a 
factor of 3. (Individual wave amplitudes are considerably larger.) Fig. 22a shows at 60 km that the 
oscillations extend from the equator to high latitudes and occupy primarily the winter hemisphere. At 
90 km (Fig. 22b), however, the pattern is markedly different. The oscillations populate the region 
around the equator and extend at higher latitudes into the summer hemisphere. At both altitudes, the 
amplitudes approach values close to 30 m/s. 

Our analysis indicates that the changes in the seasonal pattern of the m = 3 PW, also seen in m = 4, 
are produced by variations in the baroclinic instabilities that prevail in the stratosphere and 
mesosphere. This is shown in Fig. 23, where we present at 50 and 90 km the normalized differences 
between the relative latitudinal pressure and density gradients for the zonal mean (m = 0). The 
gradients are largest at 50 km in winter (Fig. 23a), but at 90 km in summer (Fig. 23b). In the 
mesosphere, the reversal of the latitudinal temperature variations (Fig. 20) is accompanied by changing 
baroclinicity, which in turn produces the changes in the seasonal PW pattern in Fig. 22. 

6. Summary 

We summarize numerical results from the Numerical Spectral Model (NSM). The model is 
integrated from the ground into the thermosphere, and it is fully non-linear. The spectral formulation 
allows us to analyze the individual, and interacting, dynamical components with different zonal wave 
numbers, which describe the zonal-mean variations, tides, and planetary waves. 

The NSM employs the Doppler Spread Parameterization (DSP) for small-scale gravity waves 
(GW) developed by Hines (1997a, b). This formulation reproduces a broad range of dynamical 
features observed in the zonal-mean (HMCF) - in particular the quasi-biennial oscillation (QBO) of 
the stratospheric circulation, which is generated with GW momentum deposition owing to critical level 
absorption. 

The kind of GW source that generates the QBO can also amplify the tides (Walterscheid, 19861). 
This is shown in the present paper with numerical results, which demonstrate that the GW amplify the 
tides and planetary waves, and presumably the planetary-scale inertia gravity waves and related pseudo 
tides. 
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In Hines' GW parameterization, the wave momentum flux is conserved, which has important 
consequences. The GW momentum source that generates the QBO is not available to amplify the tides 
at higher altitudes. Through GW filtering, the QBO signature is transferred to the diurnal tide to 
generate inter-annual variations in the upper mesosphere. The GW momentum, expended in the 
annual variations of the mean zonal winds, produces equinoctial amplitude maxima in the diurnal tide. 
And GW interactions with the diurnal tide significantly affect the seasonal variations of the semi- 
diurnal tide. Numerical experiments demonstrate that GW produce nonlinear interactions between 
migrating tides and planetary waves, to generate non-migrating tides. 

As the GW propagate up through the atmosphere, they cairy the accumulating imprints of 
interactions with the zonal-mean variations, tides, and planetary waves. The accompanying nonlinear 
interactions, due to GW filtering and momentum deposition, produces a climatology that becomes 
increasingly more complex at higher altitudes — apparent in our model results and in qualitative 
agreement with observations. 
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Figure Captions: 


Figure 1 . Equatorial oscillations (QBO and SAO) generated with 2D (a) and 3D (b) versions of the 
NSM, applying identical GW sources. The QBO in the stratosphere (30 km), and the SAO in the 
mesosphere (55 and 80 km) have larger amplitudes in 2D. The GW amplify in 3D the tides and 
planetary waves at the expense of the QBO and SAO. (Figure taken from Mayr et al., 1 999.) 

Figure 2. Schematic diagram illustrates the processes that generate and influence the migrating and 
non-migrating tides, planetary-scalc inertio gravity waves, and planetary waves. In the present paper, 
numerical experiments are discussed that illuminate the GW interactions enumerated. 

Figure 3. Amplitude and phase of the diurnal tide generated with (a) and without (b) GW interactions. 
The GW increase the amplitudes and decrease the vertical wavelengths. (Figure taken from Mayr et al., 
1998.) 

Figure 4, Schematic diagram illustrating the GW interaction, which causes momentum deposition in 
regions of vertical wind shears, analogous to critical level absorption that produces the QBO. (Figure 
taken from Mayr et al., 1999.) 

Figure 5. (a) Temporal variation of diurnal tide at 95 km and 18° latitude shows amplitude maxima 
near equinox, (b) The mean zonal winds change direction between 50 km and 80 km due to GW 
interaction, and their seasonal variations produce nodes close to equinox. During spring and fall, the 
GW propagate up without being much absorbed, and this produces semi-annual amplitude maxima in 
the diurnal tide (a). (Figure taken from Mayr et al., 1999.) 

Figure 6. Amplitude and phase of the diurnal tide computed with and without GW interaction. (Figure 
taken from Mayr et al., 2001 .) 

Figure 7. (a) The planetary waves (PW) for m = 3 at 80 km peak around solstice, and they are 
amplified by GW. (b) Diurnal tide (DT) at 95 km computed with and without the m = 3 and 4 PW. 

The GW momentum expended on the PW is not available to amplify the DT, which produces lower 
amplitudes in summer. (Figures taken from Mayr et al., 2001 .) 

Figure 8. (a) Semidiurnal tide computed with and without GW. Like the diurnal tide (Fig. 6), the GW 
source causes amplification, (b) Semidiurnal tide computed with and without the excitation of the 
diurnal tide. During the spring and fall seasons, the GW amplify the diurnal tide at the expense of the 
semi-diurnal tide. (Figures taken from Mayr et al., 2001.) 

Figure 9. Mean zonal winds in the lower stratosphere near the equator represent the QBO with a 
period of about 24 months (a). In the zonal winds of the diurnal tide at 95 km near the equator (b), the 
QBO periodicity is prevalent (c), and it appears also at 18° latitude where the tide peaks (d). (Figures 
taken from Mayr and Mengel, 2005.) 

Figure 10. Normalized meridional (a) and vertical (b) advection terms generated by the QBO 
interacting with the diurnal tide. Pronounced QBO signatures are apparent but only at altitudes below 
about 60 km. (Figure taken from Mayr and Mengel, 2005.) 

Figure 11. (a) Running average over 364 days for the mean zonal winds at 4° latitude, (b) Normalized 
GW momentum source, analogous to Fig. 10, shows inter-annual variations with a period close to 24 
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months extending into the upper mesosphere, (c) Amplitude of diurnal tide varies in phase with the 
GW momentum source. (Figure taken from Mayr and Mengel, 2005.) 

Figure 12. Meridional winds at 100 km for the migrating and non-migrating diurnal tides, (a) The 
westward propagating tide peaks near 20° latitude, and the maxima near equinox exceed 90 m/s. (b) 
The eastward propagating tide for m = 1 is highly variable and peaks at the equator with amplitudes 
that approach 30 m/s. (c) The variable stationary tide for m = 0 is also confined to low latitudes, and 
its amplitudes can reach 40 m/s. (Figures taken from Mayr et al., 2005a.) 

Figure 13. Non-migrating tides are produced in the NSM with planetary waves (PW) that are 
generated internally by instabilities in the zonal-mean, (a) The stationary diurnal tide for m = 0 is 
produced with m = 1 PW, and correlations are evident, (b) Without zonal-mean heat source, PW and 
non-migrating tide are not generated. (Figures taken from Mayr et al., 2003a.) 

Figure 14, Illustration of proposed GW interaction for generating nonlinear coupling between 
migrating tides and PW. In the positive phase of the PW, for example, the eastward (positive) 
propagating GW are preferentially absorbed. The GW are then less effective in amplifying the 
eastward phase of the tide, which produces a non-migrating tide that is modulated by PW. (Figure 
taken from Mayr et al., 2005b.) 

Figure 15. Non-migrating m = 0 tide generated with and without GW interaction. (Figure taken from 
Mayr et al., 2005b.) 

Figure 16. Computer experiment with standard and time independent m = 0 GW source, (a) The 
migrating tide (m = 1) is not affected, (b) With constant GW source (lower panel), the PW amplitudes 
are larger, (c) Non-migrating tide generated with migrating tide (a) and PW (b). The time dependence 
of the GW source produces a much larger tide, although the PW involved are somewhat weaker. This 
demonstrates the importance of GW for the non-linear coupling between migrating tide and PW. 
(Figures taken from Mayr et al., 2005b.) 

Figure 17. Oscillations for m = 1 at high latitudes are shown with periods around 1 0 hours. Like tides, 
the waves grow with height and have vertical wavelengths of about 20 km. The amplitudes of 
meridional winds (a, b) and temperature variations approach in the upper mesosphere 40 m/s and 5 K, 
respectively. (Figures taken from Mayr et al., 2003b, 2004b.) 

Figure 18. Meridional wind spectrum for m = lat 100 km shows westward propagating 12-hour semi- 
diurnal tide. Oscillations around 9 and 10 hours propagating westward at polar latitudes, and 9-hour 
waves propagating eastward at mid latitudes. (Figure taken from Mayr et al., 2004b.) 

Figure 19. (a) Seasonal variations of westward and eastward propagating 10-hour waves show peak 
amplitudes during winter and spring seasons, (b) Without excitation of diurnal tides, 10-hour waves 
are generated with comparable amplitudes, but the seasonal variations are different. (Figures taken 
from Mayr et al., 2004b.) 

Figure 20. (a) Temperature variations for December solstice show a pattern in qualitative agreement 
with observations, (b) Reversing latitudinal temperature variations, and related meridional circulation, 
are conducive for generating baroclinic instabilities in the troposphere and middle atmosphere, which 
produces the planetary waves in the model. (Figure taken from Mayr et al., 2004a.) 
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Figure 2 1 . Planetary waves at 90 km and mid latitudes for m = 1 and 3, generated with tropospheric 
source (a) and without (b). (Figures taken from Mayr et al., 2004a.) 

Figure 22. Planetary waves for m = 3 at 60 and 90 km show changing patterns of seasonal variations. 
Peaks occur in winter at 60 km, and in summer at 90 km. (Figures taken from Mayr et al., 2004a.) 

Figure 23. Meridional gradients of m = 0 pressure and density are measures of baroclinic instability. 
Peaks occur in winter at 50 km (a) but in summer at 90 km (b), which is consistent with the changing 
planetary wave patterns in Fig. 22. (Figures taken from Mayr et al., 2004a.) 
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Gravity Waves (GW) 


(a) GW interaction produces temperature and 
wind reversals that lead to instabilities 

(b) GW filtering produces seasonal variations 

(c) GW interaction amplifies tides 

(d) GW source of QBO 


(e) GW filtering produces inter-annual variations 

(f) GW source amplifies planetary-scale inertio gravity waves 

(g) GW' interaction amplifies planetary waves 

(h) GW filtering contributes to nonlinear interaction between 
migrating tides and planetary waves 
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Gravity Wave Interaction with Tide 
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Diurnal Tide: m = 0, Meridional Wind, 100 km, 18° N, 3 day window 
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